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Abstract 

We continue our investigation of the configuration space of general rel- 
ativity begun in I (J. Guven and N. O Murchadha, ^r-qc/9411009| ). 
Here we examine the Hamiltonian constraint when the spatial geome- 
try is momentarily static (MS). We begin with a heuristic description 
of the presence of apparent horizons and singularities. A peculiarity of 
MS configurations is that not only do they satisfy the positive quasi- 
local mass (QLM) theorem, they also satisfy its converse: the QLM is 
positive everywhere, if and only if the (non-trivial) spatial geometry is 
non-singular. We derive an analytical expression for the spatial metric 
in the neighborhood of a generic singularity. The corresponding curva- 
ture singularity shows up in the traceless component of the Ricci ten- 
sor. As a consequence of the converse, if the energy density of matter is 
monotonically decreasing, the geometry cannot be singular. A super- 
metric on the configuration space which distinguishes between singular 
geometries and non-singular ones is constructed explicitly. Global nec- 
essary and sufficient criteria for the formation of trapped surfaces and 
singularities are framed in terms of inequalities which relate some ap- 
propriate measure of the material energy content on a given support 
to a measure of its volume. The sufficiency criteria are cast in the 
form: if the material energy exceeds some universal constant times 
the proper radius, £q, of the distribution, the geometry will possess 
an apparent horizon for one constant and a singularity for some other 
larger constant. A more appropriate measure of the material energy 
for casting the necessary criteria is the maximum value of the energy 
density of matter, pMax^ if Pmax^o < some constant the distribution 
of matter does will not possess a singularity for one constant and an 
apparent horizon for some other smaller constant. These inequalities 
provide an approximate characterization of the singular/non-singular 
and trapped/non-trapped partitions on the configuration space. Their 
strength is gauged by exploiting the exactly solvable piece-wise con- 
stant density star as a template. Finally, we provide a more transpar- 
ent derivation of the lower bound on the binding energy conjectured 



by Arnowitt, Deser and Misner and proven by Bizon, Malec and O 
Murchadha and speculate on possible improvements. 



1. INTRODUCTION 

This is the second paper in a series in which we examine the solution of the constraints 
in general relativity under the restriction that both the spatial geometry and the extrinsic 
curvature be spherically symmetric. In the first paper (hereafter referred to as paper I), 
we examined various universal features of the constraints [1] . In particular, it was shown 
following [2] how the constraints can be recast in a symmetrical form with respect to 
the optical scalars of the theory. We introduced an appropriate large class of foliations 
which are in some sense the natural foliations of a spherically symmetrical spacetime. We 
discussed the positivity of a quasi-local mass and related local bounds on the canonical 
variables of the theory. Finally, a tentative realization of the configuration space was 
proposed. 

In this paper we will specialize to solutions of the constraints when the geometry is 
momentarily static (MS) and asymptotically fiat.* The extrinsic curvature and the radial 
material current now both vanish, Kab = = J and the radial momentum constraint 
becomes vacuous. We follow the notation established in paper I. The proper (or geodesic) 
radius £ is treated as the independent radial coordinate in our parametrization of the spatial 
geometry. The Hamiltonian constraint then provides an ODE for the circumferencial radius 
as a function of £. We defer the examination of non- vanishing extrinsic curvature to a third 
paper (III) [3]. 

Whether the solution remains static, subject to the initial conditions Kab = so that 
Kab = will depend on the stresses which only show up in the dynamical equations. What 
is true is that if one is to construct globally static solutions, they had better satisfy the 
Hamiltonian constraint with Kab = to begin with. 

While the constraints simplify enormously when Kab = 0, many if not all of the sub- 
tleties involved in the solution of the initial value problem obligingly manifest themselves. 
This makes such configurations important as a theoretical proving ground. 

Another feature of these solutions which makes them interesting is that they are the 
only geometries which simultaneously satisfy the Lorentzian and Euclidean constraints. 
In the semi-classical approximation to tunnelling, they correspond to the classical turn- 
ing points. They therefore represent the boundaries between the classically allowed and 



* A static solution is a solution possessing a surface forming timelike Killing vector. It 
is easy to show that the spacelike surface orthogonal to the Killing vector satisfies Kab = 0. 



forbidden regimes [4]. 

Our goal will be to provide a characterization of the configuration space, C, of all spa- 
tial three-geometries satisfying the Hamiltonian constraint essentially by whatever means 
we can. From one point of view, of course, C could be characterized by whatever material 
degrees of freedom we introduce and, as such, is trivial even if the corresponding Hamil- 
tonian is not [5]. If matter is modelled by a scalar field, for example, C will simply be 
all spherically symmetric scalar fields with an appropriate falloff. What is ignored here, 
however, is that to each point in the configuration space there also corresponds a spa- 
tial geometry. In particular, some configurations will correspond to geometries possessing 
apparent horizons, with a proper subset of these corresponding to geometries possessing 
singularities. The boundaries in C partitioning geometries with apparent horizons or sin- 
gularities from those which do not are important landmarks on this space. 

In Sect. 2 we discuss generic analytical properties of the Hamiltonian constraint. We 
choose, as our independent data, a spherically symmetric energy density distribution. This 
consists of a non-negative function p, given as a function of the proper radius £ on the 
positive real axis, [0, oo). The Hamiltonian constraint is then solved to give the spatial 
geometry. This is to be regarded as a kinematical quantity, entirely determined by the 
source distribution. We are not fussy about the differentiability of p, we would be happy 
with C^ where A^ is any integer in the interval [0, oo] . 

We suppose that the asymptotically fiat geometry is regular at its center. Regular 
solutions of the Hamiltonian constraint are determined uniquely once p is specified. For 
some choices of p, however, we find that the circumferential radius R goes to zero at a 
finite value of £ = is- 

Let us suppose we scale some small p, say po. Initially, the circumferential radius 
R increases monotonically with the proper radius. As the scale is raised a critical point 
is always reached in which the geometry develops an apparent horizon. When Kab = 0, 
this corresponds to the development of an extremal embedded two-sphere. At this critical 
scaling, the horizon will be degenerate. As p is scaled further, the horizon will bifurcate 
into an outer minimal two-sphere and an inner maximal two-sphere. The interpolating 
two-spheres which foliate the spatial region between the two will all be trapped. Typically, 
if Po is scaled further, the geometry will become singular at some other critical scaling 
corresponding to the pinching off of the minimum two-sphere to i? = at £ = £s- AH points 
on this two-sphere would then be identified, the geometry would possess a bag of gold at its 



center. This is the only way that the spatial geometry could possess a singularity. These 
three-dimensional singularities are the only ones we discuss. These form only a subset 
of the possible four-dimensional singularities that could arise when regular initial data is 
evolved using the dynamical equations to find a solution of the Einstein equations. One 
particular class of singularities (which are simultaneously three- and four-dimensional) we 
do not discuss are those which arise when p at a point becomes unboundedly large. 

Physically, it might be argued that the singular geometries we consider are irrelevant 
from the point of view of the asymptotically flat spacetime as they are disconnected from 
it. It could also be argued, however, that this feature is an artifact of the spherically 
symmetric model. If the symmetry is relaxed, the singular structures which arise will, 
typically, not disconnect the spatial geometry. This is not, however, our justification for 
considering singular structures. They are important for the reason indicated by the scaling 
argument. An infinesimal change in p can convert a regular geometry into a singular one. 
The description of the configuration space is incomplete if such singular geometries are 
ignored. 

A remarkable feature of spherically symmetrical geometries is the positivity of a well 
defined quasi-local mass (QLM) whenever the spatial geometry is regular and the sources 
satisfy the dominant energy condition [1,6]. Its converse is not, however, generally true — 
there are singular solutions of the constraints with a positive everywhere quasilocal mass 
[3] . A peculiarity of MS configurations is that they do satisfy the converse of the positive 
QLM theorem. If the QLM is positive everywhere, the geometry is necessarily non-singular. 
We can now exploit the definition of the QLM to derive an analytical expression for the 
local form of the spatial metric in the neighborhood of a generic singularity. Surprisingly, 
the QLM remains finite at the singularity. The constraints guarantee that the scalar 
curvature is also finite everywhere, even if the geometry is singular. If there is a curvature 
singularity, it must show up as a divergence in the traceless component of the Ricci tensor. 
Non-generic radial density profiles exist in which the geometry degenerates without any 
curvature singularity where R pinches off. The interior region behind the singularity can 
be interpreted as a regular closed universe. 

Our understanding of the analytic form of the metric in the neighborhood of a sin- 
gularity, facilitates the identification of a supermetric on C in the manner of DeWitt [7]. 
Specifically, in Sect. 3 we identify a non-ultralocal metric which assigns a finite norm to reg- 
ular geometries and an infinite norm to singular geometries. With respect to this metric, 
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the boundary in C separating singular geometries from non-singular ones can be consigned 
to infinity. 

A surprising consequence of the converse of the positive QLM theorem is that if 
p is strictly monotonically decreasing, the geometry cannot be singular (if it is regular 
at the origin). What this demonstrates is that p can be scaled by an arbitrarily large 
factor on a given support and yet the geometry remain regular everywhere. This counter- 
example demonstrates clearly that the scaling argument breaks down as a description of 
the development of singularities. For while an apparent horizon always develops at some 
critical scaling, a singularity does not necessarily develop. Such geometries may not be of 
great physical interest because we would expect that any classical singularity that might 
form in such a collapsing configuration would form at its center where we have assumed 
that the geometry is regular. However, monotonicity does provide a useful boundary in C 
separating geometries with very different properties with regard to singularity formation. 

We will illustrate these points explicitly in Sects. 4 and 5 using a piece-wise constant 
density star. 

We discuss the constant density star in greatest detail. The model does not admit 
curvature singularities; its only pathologies are metric degeneracies. Because of this one 
needs to be extra wary of making generalizations. We illustrate how a mechanical analogue 
can be used to study the solution space. 

The local measures of energy and size which are most appropriate for measurements 
performed outside the star consist respectively of the QLM m and R at the surface of 
the star. What is clear already in this model, however, is that m provides a very poor 
measure of the material energy content just as R becomes a poor measure of the size when 
the geometry is highly curved. In particular, m (or R ) may be infinitesimally small in 
geometries possessing apparent horizons. 

It is true that an observer living outside the star does not have any choice but to 
exploit these 'local' measures. Our goal, however, is to chart the configuration space, 
above all, in the neighborhood of configurations with apparent horizons and singularities. 
Because we are not asking the kind of question a remote orbiting observer in one of these 
geometries might ask, there is no reason why we should handicap ourselves with variables 
more appropriate to another scenario. The geometry behind a horizon is not out of bounds. 
In is in this spirit that we identify our global measures of energy content and size. 

A measure which does appear to correspond quantitatively to the material energy 



content is the integrated energy density over the spatial volume, M. The most appropriate 
measure of size is the proper radius of the support of matter, Iq. The monotonicity of these 
two variables means that they fare better than m and R in the strong field region we are 
interested in. We examine in detail the dependence of M on £o- 

A two-density model is discussed in Sect. 5. We demonstrate explicitly that when the 
inner density exceeds the outer one, the model can never possess a singularity no matter 
how great the density or large its support. If, on the other hand, the outer density is the 
greater, singularities can appear if the star is large enough. If the outer density is tuned 
appropriately, the curvature singularity becomes a metric degeneracy. 

In Sect. 6 we follow Bizon, Malec and O Murchadha (BMOM) by concentrating on 
the formulation of global necessary and sufficient conditions determining when the spatial 
geometry will possess apparent horizons and singularities [8,9,10]. The motivation behind 
this work was Thome's hoop conjecture [11]. These conditions are cast as inequalities 
involving appropriate global measures of the energy content of the matter distribution and 
of its support. 

The sufficient conditions can be cast in the form: if M > C Iq, where C is a constant, 
then the star possesses a trapped surface for one constant, and a singularity for some 
larger constant within the sphere of proper radius Iq. On the subset which corresponds to 
a monotonicically decreasing density these inequalities can be strengthened. However, the 
latter inequality now provides a universal bound on M. To understand how this comes 
about suppose we are given some monotonically decreasing density, po- While one might 
imagine that by scaling po by some large constant we could make M arbitrarily large, 
this is misleading because the volume element appearing in M contains R^. Given po? 
the constraint determines the value of R. When we solve the constraint with the scaled 
value of p, the net effect is to decrease R. This can occur in such a way that M, in fact, 
saturates. 

In the original work of BMOM, conformally ffat coordinates were exploited. The 
attraction of these coordinates is that they do not rely on the spherical symmetry of 
the problem. The super-Hamiltonian constraint reduces to a non-linear elliptic PDE on 
ffat R^. This is offset, however, by the fact that the background geometry is no longer 
the physical geometry and one must contend with unphysical bifurcations in the solution 
space of the elliptic equation. What is more, the existence of the inequalities discovered 
by BMOM is not obvious when the constraints are expressed in this way. Their derivation 



demonstrated a considerable level of ingenuity. By expressing the constraints directly in 
terms of physical quantities, we are able to provide very simple alternative rederivations 
of the BMOM inequalities which are both more economical and transparent than theirs. 

The shortcomings of M become apparent (even in the constant density model) when 
we attempt to provide necessary conditions on the specification of the boundaries parti- 
tioning the configuration space. For, li M < C Iq, C a, constant, on the surface of the star, 
for some constant no matter how small, even if it is true that the surface is neither trapped 
nor singular, this is not true of the interior. What is worse, we find it impossible to provide 
even this weak form of necessary condition involving the surface geometry when currents 
fiow [3]. This contrasts dramatically with the relatively straightforward generalization of 
the sufficiency conditions in this context. 

A more appropriate measure of the energy content of the star for phrasing the neces- 
sary condition is provided by the maximum value of the energy density, pMax- In Sect. 7, we 
demonstrate that if pMax^o < constant, then the interior is free of trapped surfaces for one 
constant, and free of singularities for some other larger constant. The proof requires the 
introduction of a battery of simple Sobolev and other inequalities. These bounds possess 
extremely non-trivial generalizations to J 7^ 0. 

As a demonstration of the effectiveness of our method, in Sect. 7 we rederive the 
BMOM proof of Arnowitt deser and Misner's conjectured inequality placing a lower bound 
on the binding energy of a compact spherically symmetric system [10]. 



2 GENERIC ANALYTIC FEATURES OF SOLUTIONS 

The Hamiltonian constraint provides a second order non-linear ODE for R which is 
regular everywhere except at R = 0* 

RR" -hl-{R'f)+4TvR^p = 0. (2.1) 

The closure of the geometry at its base point -R(O) = 0, is the only boundary condition we 
are free to impose. We will require that the geometry be regular at this point [1], so that 

R'{0) = ±1. (2.2) 

Otherwise, it is clear that R" would be infinite and the geometry singular at the origin. R 
can be either positive or negative — the physical metric only depends on R^. However, the 
only way R can change sign is by passing though a singularity or degeneracy with R = 
at some non-vanishing value of £. If the geometry is regular, the sign of R will not change. 
Our convention will be to choose the positive sign at £ = 0. 

Because Eq.(2.1) is singular at the origin, we need to differentiate the equation n — 1 
(not 71 — 2) times to determine the n^^ derivative at this point, R^'^\0). When we do this, 
we obtain the power series expansion for R{£) in the neighbourhood of £ = 0, 

i?(£) ^ i _ l![^£3 ^ (^(^4) _ (2.4) 

9 
The most convenient way to solve Eq.(2.1) is to exploit the fact that the equation 
possesses a first integral equating the quasi-local mass defined at a given value of £ in 
terms of the interior energy content. In paper I, this was done in a completely general 
gauge invariant context in which neither Kab nor J need be zero [1] . If we set both Kab 
and J equal to zero in Eqs.(4.5) and (4.6) of [1] we obtain^^ 



* The Hamiltonian constraint will also assume this form in either of the gauges, Kji = 

or Kc + 0.5Kii = which mimic a MSC (see paper I). In the latter case, however, apparent 

horizons will not coincide with minimal surfaces and we need to be careful when we carry 

over the MSC analysis. 

"1^ That this is also true in the polar gauge discussed in paper I is obvious. It is also true 

in the other mimicking gauge {a = 0.5) but we must remember that what we call 'm{£) in 

Eq.(2.5) is no longer the QLM. 

10 



W^ = l-^. (2.5) 

where m{i) is given by 

m{i) =Atx MR^R'p. (2.6) 

Eqs.(2.5) and (2.6) are a first integral of Eq.(2.1) which implement the condition of regu- 
larity at the origin, Eq.(2.2), explicitly. 

Outside the support of matter m{l) coincides with the ADM mass, rrioo and 

{R'f = 1 - ^ ' (2-50 

If p = everywhere, to = and space will be fiat, R{1) = £. If nioo is positive and R 
is large the second term on the right hand side of (2.5') becomes negligible and R ^ i 
to leading order. Asymptotic fiatness is guaranteed by the constraint. The 'boundary 
condition' at infinity does not need to be prescribed by hand. The ADM mass is encoded 
in the next to leading order, 

R^i-mooln£. (2.7) 

Because nioo is positive R' must tend to one from below. 

The power series (2.4) would appear to suggest that, once regularity at the origin is 
implemented, the solution will remain regular (regardless of the singularity of the differ- 
ential equation aX R = 0) so long as p{£) remains finite. Certainly, if p and its support are 
small, on physical grounds we would expect the solution to be regular with R increasing 
monotonically with i between £ = and infinity. If p (or its support) is sufficiently large, 
however, R{£) need not continue to increase monotonically with £. We can imagine doing 
this by scaling on some initial small value. Beyond some critical scale factor, the curvature 
of the geometry will be sufficiently large that R' will vanish at some value of £, an apparent 
horizon forms, outside which R{£) begins to decrease. Note that this behavior is consistent 
with the convexity of R{£) given by (2.4) in the neighborhood of the origin. We will be 
more precise below. 

What can go wrong beyond this point is that R might not recover but will continue 
decreasing until it reaches zero (where the constraint equation is singular) at some non- 
vanishing value of £. If we extrapolate the first two terms in Eq.(2.4) to larger £, we see 
that R{i) returns to zero when p{0)£q = 9/47r — in reasonable agreement with the value 
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for the constant density star discussed in Sect. 4. The important point is that this is the 
only way that the spatial geometry can be singular.* 

In Sect. 5 of paper I, we proved the positivity of the QLM, m: if the geometry is regular 
and the weak energy condition is satisfied everywhere then R'^ < 1 everywhere and, as a 
consequence, m > 0.^ 

Obviously, m{£) goes to zero as one approaches the origin. It does so faster than R 
so that R' approaches unity. Starting from the origin and moving outwards, we see that 
m{£) = and R' = 1 until we run into some nonzero p. If the support of p is bounded away 
from the origin we get a locally flat region around the origin. Once we run into nontrivial 
p we have that m{£) must increase to reach some positive value. If m{£) is bounded from 
below by some positive constant value for all positive values of £ outside a neighbourhood 
of the origin then R{£) cannot become small in this region if Kab = 0. This is because this 
would imply that the second term on the right hand side of Eq.(2.5) must exceed the flrst 
at some point which is impossible. Thus the geometry cannot become singular. 

We can therefore claim that 'm{£) > everywhere except at the origin or in a neigh- 
bourhood of it if and only if the geometry is non-singular. This neighbourhood, of course, 
covers the whole space only in the trivial case where p = and we have flat space. The 
set of geometries with m{£) = somewhere away from the origin is a set of zero measure 
and will be treated separately below after we have discussed the constant density star. 

What is the analytical description of a generic singularity? Suppose that R'{£) < —1, 
or equivalently m{£) < at some point £. Eq.(2.1) tells us that R" < 0. Hence R' must 
become more negative. R must return to zero at some flnite value £s beyond £, R{£s) = 0. 
If p is bounded on the interval [0,^^] it is clear that m{£s) reaches a flnite negative value. 
Let us assume that R goes to zero like (£5 — £)'>', where 7 is some positive constant. We then 
expect R' to behave like {£s — £)^~'^- From Eq.(2.5) we see that m{£) goes like {£s — £)^^~'^ 
and that m'{£) should go like {£3 — £)^''~^. However, if we differentiate Eq.(2.6) we flnd 
that m'{£) equals AtiR'^R' p and this means it should behave like {£3 — £Y'^~^. This is 
inconsistent. The only way out is if 87 = 2. This means that m{£) goes to a constant and 



* Even when p is large, R{£) cannot diverge at any flnite value of £. This is because this 

would require R' to diverge whereas the only place it can diverge is at the singular point 

of the differential equation i? = 0. 

t The contrapositive of this statement is that if m{£) < at some £, the geometry must 

become singular at some point beyond £. 
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when we differentiate we do not get the 87 — 3 term. 

Near £5, the second term in Eq.(2.5) will dominate so that 



(i^')^^-^. (2. 



This equation can be integrated to give 



3,2/3 



R^l-] (-2m{is)f'{is-if/', (2.9) 



entirely consistent with the argument above. R' becomes unbounded at is as — (£5— £)~^/^. 
If, therefore, R'^ > 1 {R' < —1) anywhere it must subsequently diverge. It is clear, from 
(2.6) that not only does m remains finite at £5, so also does the integrand pR^R' . 

All points with £ = £s are zero distance apart so they must be identified. The 
geometry now contains one of Wheeler's bags of gold [12] behind the singularity. Because 
the singularity cuts off the interior geometry, it is customary in the present context to 
discard the solution as unphysical. The exterior does, however, also contain this singularity. 

Generically R' will diverge as we approach a 'singularity'. We will call this a strong 
singularity. If, however, we fine tune the density so that m{£s) = then R'{£s) = —l? 
and R"{£o) = 0. The bag closes smoothly, and there are no singularities in the curvature. 
We can always approximate such a geometry arbitrarily closely by a regular geometry. We 
will refer to this as a weak singularity. 

A beautiful consequence of the converse of the positive QLM theorem is that it permits 
us to identify an important class of non-singular geometries. We note that, by a simple 
integration by parts, 'm{£) can be expressed in the alternative form 

47r r^ 

m{£) = —R^p-An d£R^p'. (2.6') 

If the energy density is monotonically decreasing, p' < 0, then 

Att 
m{£)>—R^p. (2.10) 

o 

In particular, if p is a monotonically decreasing function of £, and it is strictly monotonically 
decreasing in some neighbourhood of the origin then m{£) > 0, and so the geometry is 
necessarily non-singular. 
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We note that neither the converse of the positivity of the QLM nor the consequence 
survive when J 7^ 0. In general, a positive value of m{£) everywhere does not guarantee 
the non-singularity of the geometry. 

What is the geometrical nature of a generic spatial singularity? Even though the 
geometry might be singular, we assume p is always finite. The constraint therefore guar- 
antees that the scalar curvature remains finite. The metric singularity must therefore show 
up in the Ricci tensor. We can express 

'Tlab = '^/:nanb + TlR{gab-nani,), (2.11) 

in the same way as we did Kab- AH curvature scalar s can be expressed in terms of 7^£ and 
TZji. In appendix I, we show that 

i T (2.12a, 6) 

While the scalars Re and Rr both diverge as one approaches a strong singularity, the 
sum IZ = IZc + 27?.jj remains finite. In fact, as we approach a strong singularity, Tic,R ~ 
±i?'^/i?^ ~ ±(^s — £)~^. The singularity need not occur on the support of matter. At a 
weak singularity IZab remains finite. 



3 TO THE SINGULAR BOUNDARY OF CONFIGURATION SPACE 

3.1 Supermetrics 

In this section we construct a supermetric on the configuration space, C, of spheri- 
cally symmetric asymptotically fiat geometries satisfying the Hamiltonian constraint at a 
moment of time symmetry. Such geometries are completely characterized by the function, 
p{t). We can introduce a fiat line element on C 

2 / jfi s „l f)\2 



\\bp\\'=\ dibpity. (3.1) 

which is ultra-local with respect to p. The fiatness of the metric makes life very convenient 
in that it is then simple to determine the 'distance' between finitely separated geometries. 
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The shortcoming of the metric, (3.1), is that it does not tell us anything about the structure 
of the underlying geometry. 

To do this, we need to construct a line element which depends explicitly on the 
geometry. In the full theory, the commonly accepted choice is the ultra-local DeWitt 
supermetric 



ll'^^llLvKo = Jd'x^G'^'^''6gaBSgcd. (3.2) 



where G"'^ '^^ = | ^gacgbd _|_ gadgbc^ _|_ Qgahgcd ^^^ Q jg gome constant which is determined 
by the requirement that G"'^ ^'^ be positive definite. On the subset of spherically symmetric 
geometries, ||5-R||DeW0 is given, up to an irrelevant constant, by 

nOO 

\\SR\\leWo= di5R{i)\ (3.3) 

Jo 

which, remarkably, is also fiat. Note that the DeWitt supermetric is not itself fiat. The 

metric it induces on the subspace of spherically symmetric geometries is flat. 

We note that with respect to any other parametrization of the spatial metric, the form 

assumed by the DeWitt supermetric will be more complicated. For example, with respect 

to the Schwarzschild parametrization (when valid) in the notation of I: 

WomOeWO- dK 2 ■ ('^•4J 

Unfortunately, no matter how we parameterize it, the norm on any asymptotically geome- 
try, even the flat geometry, is inflnite with respect to the ultra- local metric Eq.(3.2) simply 
because i? ~ £ at inflnity. In fact the distance between any two geometries with differing 
values of TO will also be inflnite with respect to ||5i?||DeW- The reason for this is that the 
asymptotic behavior given by Eq.(2.7) implies that i?i — i?2 ~ (?tIoo1 — "^002) ln£ which 
is not square integrable. It is not, however, clear if this is anything to worry about, moo, 
after all, is conserved and, classically at least, different values of m correspond to subsets 
of C which do not intersect. Even quantum mechanically, it would appear that the theory 
must be constructed with a flxed value of moo • 

It is, however, simple to construct a metric which gives a flnite norm between ge- 
ometries with different moo's. We do this by relaxing the criterion of non-ultralocality 
by considering a metric involving a flnite number of derivatives of R. If this seems like a 

15 



high price to pay, one should note that in terms of p, even the DeWitt ultralocal metric is 
extremely non-local. Consider the metric 

/•OO 

\\SR\\leWi= / di\5R'{e)\\ (3.5) 

Now R'l — i?2 ~ ("^ooi — "^002)/^ which is square integrable. Even still, the norm of any 
geometry is still infinite with respect to the non ultra-local metric we have introduced. 
Furthermore, none of the metrics we have considered so far, discriminates against the kind 
of singular geometries which show up as solutions of the constraints. In particular, it is 
clear from Eq.(2.9) that R' is square integrable over any compact interval, whether the 
geometry is singular or not. The simplest metric that solves both problems is one which 
is second order in derivatives * 

/■OO 

\m\\leW2= di\SR"{£)\\ (3.6) 

Now ||5-R||DeW2 is finite on any non-singular geometry that satisfies the constraints. 
Eq.(2.9) indicates that, with respect to this metric, not only are strongly singular geome- 
tries pushed out to infinity, they are also rendered infinitely distant from any non-singular 
geometry. What is also very attractive is that, with respect to this metric, the norm of 
the fiat geometry R{i) = £ is zero. However, this might be disconcerting for those who 
advocate a zero geometry as the ground state. While the fiat metric in p also does this, it 
fails to penalize singularities. 

It is possible to introduce a metric which is geometrically more satisfying, which is also 
second order in derivatives by casting it directly in terms of differences in the curvature 
scalars. Unfortunately, the price we pay is that the resulting metric is no longer fiat. In 
addition, the simplest such choice, 

<>oo 

||57^|p= / diR^\6n{i)\'^ , (3.7) 

Jo 

is a step back because it fails to discriminate against singular geometries. This is solved 
by considering the metric (C = in the DeWitt supermetric) 



* In general, on any compact support, these metrics satisfy the sequence of (Poincare) 
inequalities: 



||'^-R||DeiyO < ||'^-R||DeWl < | I'^-R] IdcVK 2 < 
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/■OO 

||5Ricc||2=/ MR^G''^'"^5nabSncd 

"^° (3.8) 



oo 







|57^£(£)|2 + 2|57^fl(£)p 



3.2 The Singular Boundary 



Perhaps, the most important partition of C is into geometries which possess one or 
more singularities and those which do not, Cq. Cq is clearly open. Within Co, lie all 
geometries corresponding to (strictly) monotonically decreasing p(£). With respect to any 
of the metrics we have introduced on C, two stars with monotonically decreasing energy 
density profiles can be arbitrarily far apart. Cq is therefore unbounded, at least in some 
directions, with respect to this metric. The boundary S appears then to be disconnected. 
Our task will be to begin to characterize it. 

It is not clear if we should reject singular geometries as acceptable points on the 
configuration space. Should we treat the surface separating singular geometries from non- 
singular ones as the boundary of the physical configuration space? If this is so then we 
are faced with a major problem. Technically, it appears unlikely that we will ever be in a 
position to identify S exactly. In the quantum theory, we would be required to implement 
a boundary condition on the wavefunction on S. What is to be hoped for is that the 
potential determining the dynamics in the configuration space will be sufficiently steep in 
almost all directions as we enter the grey area in the neighborhood of this boundary that 
wave-packets will get reffected back onto the support of 'good' geometries. An analogy 
which might be useful is the manner in which all wave packets which enter the Taub 
channels in homogeneous mini-superspace quantum gravity appear to get reffected back 
into the generic mixmaster mini-superspace [13]. 

There is another natural partition within Cq, with boundary T into geometries which 
possess trapped surfaces and those which do not, Vq. One would expect that Vq is bounded 
with respect to any reasonable norm. A geometry can only be singular if i?' = somewhere 
so that it also possesses a trapped surface. That there is no natural supermetric which 
discriminates against geometries possessing an apparent horizon indicates the very different 
natures of the boundaries, S and T. 
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4. THE MOMENTARILY STATIC CONSTANT DENSITY STAR 
4.1 The Exact Solution 

In general, we cannot solve Eq.(2.1) exactly for arbitrary p{£). What we can do is 
solve a few simple exactly-solved problems consisting of piece-wise constant values of p. 
Such models can clearly approximate a generic energy profile arbitrarily closely. We can 
therefore exploit them as a guide to identifying the generic features of solutions. 

The simplest model there is is the constant density star.* If p is a constant in some 
interval i < io containing the origin, we can recast Eqs.(2.5) and (2.6) as the single 
equation, 

R'^ + ^r2 ^ 1 (4.1) 

in this region. Outside the support of matter, 

R" + ^ = l, (4.2) 

where rrioo = ^P-R(^o)^ is the ADM mass. We note that both R and R' are continuous 
across £ = £q. 

We can analyse Eqs.(4.1) and (4.2) by drawing on the analogy with the one dimensional 
motion of a particle {£ is time) with unit 'energy' in a time dependent potential, V{R,£), 
which changes from a harmonic oscillator to a Coulomb potential at £ = £o- 

This analogy will be particularly useful below when we examine less trivial models. In 
the present case, the particle starts out at the origin R = at £ = with unit 'velocity' 
R' = 1. The classical turning point marking the position of an apparent horizon (if one 
exists) is given hj R' = or R^ = 3/(87rp). This represents the Schwarzschild radius of 
the inner apparent horizon. Inside matter, R = R{£) is given by 

dR „ . —A I R 



l^_(u^ ^R^ 

rtc 



i?eSin-M— , (4.4) 



* We should point out that while the concept of a constant density does not depend on 
the spatial system of coordinates it does depend on the foliation. 
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or R = Re sin ( ^ 1 . where ic is given by £c = §-Rc- This is simply the metric on a round 
three-sphere S^ of radius Re- We note that the interior of the constant density star is an 
Einstein geometry, TZab = ^7^9 ab- Thus the non-singularity of TZ guarantees that of TZab- 
This is the unique such geometry. This is because all Einstein manifolds must possess 
constant TZ in dimensions greater than two (for otherwise it would be inconsistent with 
the contracted Bianchi identities) and in (i = 3, the unique spherically symmetric solution 
with positive TZ and with the given topology is part of the round S^. 

The qualitative behavior of the solution depends on the size of the star, Iq. 
If £o < ^C5 R increases monotonically both inside and out and there is no apparent horizon. 

if £o = ^c, an apparent horizon forms on the surface of the star. R is monotonic everywhere. 
We note that R" is discontinuous here. Generically, we would expect both R' and R" to 
vanish at this bifurcation. 

If £o > ^c, R increases monotonically inside the star until it reaches a maximum at the 
(inner) apparent horizon at i = £c at which R = Re- Beyond £c, R begins to decrease 
inside the star. 

If £o > 2ic, R returns to zero at i = 21 ^ inside the star at which point R' = —1. There is a 
bag of gold singularity contained within the support of matter. This is the only way that 
the geometry can exhibit a singularity in the uniform density model. The model does not 
therefore display generic singular behavior. This is its principle shortcoming. 

If 2£c > io > ic, we exit the star before a singularity is reached. However, R' is continuous 
at the boundary so that R continues decreasing outside the star until an exterior horizon 
is reached at a value of R determined hj R' = in Eq.(4b) or R = 2moo- 

We note that the constraints only admit single bubble non-singular configurations. 
Any higher number of bubbles is necessarily singular. 

4.1.1 Trapped Surfaces and Singularities: Necessary and Sufficient Conditions 

What is most interesting about this model is that it indicates that there exist bounds 
on the product of the energy density by the square of the maximum support of a given uni- 
form distribution of matter beyond which the associated geometry (i) possesses a trapped 
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surface and (ii) is singular. 

The surface of the star first becomes trapped (an apparent horizon forms) when Iq = 
7ri?c/2 so that pi^ = ||. If p^g — |f then the star wiU contain a trapped surface (either 
the surface of the star wiU be trapped or a trapped surface is contained in its interior) and 
conversely.* 

The surface of the star will be singular when Iq = tvRc/2 or pig = ■^. If pl^ > ^ 
then the star will contain a (weak) singularity (either on its surface or in its interior) and 
conversely. We summarize, 

P^^>§ iff (i) ; P^o>f iff (ii)- (4.5) 

The necessary and sufficient conditions formulated at io not only coincide but they also 
make a non-trivial statement about the interior physics. 

4.2 M vs. io, m vs. io and the characterization of trapped Surfaces and 
singularities in terms of these variables 

Unfortunately, the unambiguous partitions of the configuration space represented by 
Eq.(4.5) are a very special feature of the constant density star. If the density is not 
constant, it will not be possible in general to phrase both the necessary and the sufficient 
condition with respect to the same set of variables. When it is possible, the inequalities 
analogous to Eq.(4.5) typically will differ. The if statement will be stronger than the only 
if statement. In addition, our ability to make non-trivial statements about the interior 
based on surface measurements will, in general, depend on the choice of variables. 

We need to identify appropriate variables to characterize the material energy contained 
within a sphere of fixed proper radius when the energy density is not constant. Typically, 
this will be some (weighted) mean value of p over the interior volume. The two measures 
we introduced in paper I were the QLM m and the material energy M. 

The simplest global measure of the material energy content of a spherically symmetric 
system is 

M{io) =A-K diR^p. (4.6) 

Jo 



* We need not worry about trapped surfaces in the region exterior to the star as they 
will always be attended by an interior one. 
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Let us examine M as a function of Iq at constant p in the constant density model. We 
have 







M^Anpl deR'' = ^ 



\ Re . /2£o 



(4.7) 



What is more relevant, however, is the ratio M(£o)/^o- This is illustrated in Fig. (4.1). 

In this model it is already evident that crucial information is lost in replacing p by 
M to characterize the energy content of the star. The function i — * M(£o)/^o is neither 
surjective nor injective — it is bounded and (because the ratio oscillates) it is not unique. 
The reason for this, as we will discuss further below, is the folding of R into the definition 
of M. 

Let us attempt to formulate necessary and sufficient conditions for the presence of 
trapped surfaces and singularities in a constant density star using M{£o) as our measure 
of material energy to see how it fares. We should not be surprised that we are unable to 
do as well as we did with p — p and £q after all completely characterize the star, whereas 
M(£o) and 4 do not. 

4.2.1 Trapped Surfaces 

An apparent horizon forms at the surface of the star when M(£o) — 3£o/4. 
Sufficiency: If 

M{io) > 34/4 (4.8) 

the star will contain a trapped surface. This condition is sharp. The surface itself need 
not be trapped. However, if it is not then the interior must be singular and consequently 
must contain an apparent horizon. So far so good. 

The tightest inequality describing necessity is considerably weaker: If 

M(£o)<-^(l--), (4.9) 

4 OTT 

the star does not contain a trapped surface. An important point is that the inequality 
(4.9) describing necessity is weaker than the one describing sufficiency. In fact, inspection 
of Fig. (4.1) indicates that this inequality is not sharp and, as such, of limited value for 
there are stars with M{io) exceeding the value on the RHS of Eq.(4.9) without any trapped 
surfaces. If it had made sense to demand that the interior be non-singular in the hypothesis, 
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then, when M{io) < 3£o/4, the star would not contain a trapped surface. The possibiUty 
of singularities lurking in the interior beyond which the ratio M{£o)/£o is decreasing is 
what really limits our ability to formulate the necessary condition in terms of M{io) and 

io. 

For trapped surfaces, we can conclude that M is a good variable for the formulation 
of a sufficiency condition but not for the formulation of a necesssary one. 

4.2.2 Singularities 

The surface of the star will be singular when the ratio of M{£o) to £o is 3/4, which 
is the same as the apparent horizon ratio. This first occurs, however, when io = f^Rc/'^ 
or p£q = ^ which is different from the apparent horizon value. The description in terms 
of p and £o clearly provides a more precise characterization. When we attempt to provide 
necessary and sufficient conditions on M(£o) and £q to describe singularities we discover 
unforeseen complications. Not only are we unable to formulate conditions which coincide, 
we find that we cannot even formulate a consistent (non-vacuous) sufficiency condition of 
the form: if M{£o) > a£o for some a, then the star will contain a singularity. What we 
possess instead is a universal bound on the ratio 

M{£o)/£o<l{l + ^) (4.10) 

which holds for all geometries which satisfy the constraints, whether they are (weakly) 
singular or not. The geometry which saturates the inequality is unique. It is the configu- 
ration packing the maximum material energy into a given proper radius. It is non-singular 
but does contain a trapped surface. 

The only necessary condition we can provide with these variables is the somewhat 
trivial one implied by the corresponding trapped surface condition. If the star contains a 
singularity, it will also contain a trapped surface. However, if this criterion was poor for 
trapped surfaces, it is even poorer as a characterization of singularities. 

The reason why we are unable to provide a sufficiency condition in terms of M{£o) 
and £o which is satisfied by an non-empty set is subtle and not only a refiection of the 
inadequacy of these variables to describe non-singular geometries. We have already noted 
that a star with a monotonically decreasing p cannot be strongly singular. If the decrease 
in p is strict, neither can it be even weakly singular. We will demonstrate in Sect. 6 that 
a sufficient condition for the presence of a strong singularity is M{£o) > 2£q. We can 
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conclude that for all geometries with p' < 0, M{£o) < 21q.* The sufficiency condition 
provides a universal bound on this subset of the configuration space. The existence of the 
bound is not simply an artifact of the constant density model. 

In this regard, we should stress another distinction between the choice of p and that of 
M(£o) as a measure of the material energy which is that with p one can distinguish between 
weakly singular and non-singular geometries whereas with M(£o) one cannot. The suffi- 
ciency condition, pl^ > Stt/S is, however, unstable. For if we subject the constant density 
to a strictly monotonically decreasing perturbation, the subset satisfying this condition 
also becomes empty. 

4.2.3 m(£o) vs. 4 

If M{(.q) fares poorly for the purpose of characterizing the interior geometry, the QLM 
?b(£o) will fare even more poorly (see Sect. 7 of paper I). To demonstrate this explicitly, let 
us compare the analytical form for M, Eq.(4.7), with that for m. We have 

m(£o) = YPi?(^o)^ = §^sin=^|^. (4.11) 

?b(£o) is bounded for a given p, m < -^ = ^J?,/{?,2^lp), assuming its maximum when a 
horizon forms at Iq, decreasing thereafter as a function of Iq as matter is deposited behind 
the outer horizon. It falls to zero with the appearance of a (weak) singularity. There 
are two non-singular configurations corresponding to each value of ?Ti(^o) — one with, the 
other without an apparent horizon. 

4.3 Binding Energy in the constant density star 

In Sect. 7 of paper I, we saw that the 'binding energy', m — M is always negative. In 
the weak field limit, we note that 

M~-^ --% + ..., (4.12a) 

while 

»m-i| + .... (4.m) 



* In Sect. 7, we will show that we can do better. 
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The leading cubic terms in Iq in both coincide. It is just the 'rest' mass, rriQ = Anpl^/Z 
associated with a uniform spherically symmetric distribution of matter in special relativity. 
Thus to leading order, the difference is given by 

m-Mr^-^-^ (4.13) 

5 to 

which is the Newtonian binding energy of a spherical uniform density star. We note that 

at the other extreme where the geometry is about to turn singular, £q ~ tt-Rc, fn = and 

M ~ |M^/£o- We note that in any non-singular geometry the binding energy is bounded 

from below as follows [10] (see Sect. 7 below), 

M-m>-—. 

For a constant density star, there is the tighter lower bound, M — m > |M^/£o- The 
coefficient on the RHS increases monotonically with Iq. We conjecture that whenever 
p' < this inequality will hold. 
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5. MSCs OF THE 2-PIECEWISE CONSTANT DENSITY STAR 

The constant density model lies on the boundary of the subset of montonically de- 
creasing p, admitting weakly but not strongly singular geometries. Such weak singularities 
occur with zero measure in the set of all possible singular geometries. The constancy of p 
provides very strong global information on the surface concerning the interior. 

While it is a very useful model to get us started, it is important to recognize these 
limitations in order to avoid making erroneous extrapolations. These limitations will be 
highlighted further by examining a slighly more complicated model. 

A concrete model which does display generic behavior is a star consisting of an inner 
spherical ball of radius £i with energy density pi and an outer shell of thickness £2 — ^1 
with energy density P2: 

(pi < £ < £1 
PW = \p2 h<t<l2 (5.1) 

1 li<i- 
In the region I <i\., the solution is just one of the solutions we obtained in our examination 
of the constant density star. In £1 < £ < £2, 



where 



R' + V{R) = 1 , (5.2) 



ViR)-'-^ + '-^R\ (5.3) 



Am=^RUpi-p2), (5.4) 



i?i = i?eisin(-^j , (5.5) 

and Rli = 3/(87rpi). Am is a measure of the energy density excess in the interior. 

As we will see, the qualitative behavior of the solution in the region ii < i < £2 
depends very crucially on the sign of Am. Let us therefore examine the two signs one at 
a time (Am = is trivial). We will suppose that the interior geometry is regular in both 
cases. 

5.1 Am > 

25 



Let us first suppose that Am > 0. When Arn > 0, the potential V{R) is bounded 
from below with a positive minimum at Rm given by 

3Am^ ^/^ 



Rn.= (- . (5.1.1) 

V{R) = 1 possesses two positive roots and V{Rm) < 1 whenever Am > 0. This confirms 
that there is no critical positive value of Am beyond which no solution exists. One way to 
show this is to note that V{Rm) = 3-^ so that V{Rm) > 1 impUes Am^ > -^ — . If we 
define x = —, the inequaUty V{Rc) > 1 can be expressed in the dimensionless form 

Now, if Am > then a: < 1 which is inconsistent with Eq.(5.1.2). 

The most important feature displayed by this solution is that the geometry cannot be 
singular on the interval ii < £ < £2 or outside the star. The 'potential' provides a barrier 
keeping the solution away from R = 0. No matter how large the value of £2, the solution 
is periodic and non-singular on the interval £1 < £ < £2 oscillating between turning points 
given by the two positive roots of V{R) = 1. This provides a useful model of a (strictly 
somewhere) monotonically decreasing energy density. This should be contrasted with the 
uniform density star which turns singular beyond some critical value of £o- The inner 
higher density region has a stabilizing effect on the geometry of the outer region. 

We note that the region £ < £1 can be viewed as a momentarily static slice through 
de Sitter space. The geometry in the interval £1 < £ < £2 is the momentarily static slice 
through a Schwarzschild-de Sitter spacetime. If pi and p2 are both constant in time, 
this initial data will generate these spacetimes in their respective intervals. There will 
be, in fact, an infinite number of such geometries with the same ADM mass which we 
can label by the number of maxima of R. The existence of such geometries might have 
implications in semi-classical quantum gravity. The semi-classical calculation of tunneling 
amplitudes involves a summation over exp— 5", where S is the classical action evaluated 
on the Euclidean interpolation from one MSC to one of the infinite number of MSCs with 
the same ADM mass. This sum could potentially result in a non-trivial amplification of 
the transition amplitude. 

There is a theorem due to Schoen and Yau stating that if the volume (defined in 
terms of a certain embedded torus) supporting a given strictly positive p is made large 
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in all directions, the geometry must be singular [14]. It would appear that our ability to 
make £2 arbitrarily large without incurring any singularity is inconsistent with this result. 
However, a closer examination shows that this is not so. This is because even though the 
proper volume associated with a large value of I is large (increasing linearly with £), it is 
not large in all directions. It is essentially a cylinder of bounded radius. Thus, a large 
torus in the sense of Schoen and Yau does not fit. 

5.2 Am < 

What is more dramatic is the geometry which corresponds to Am < 0. The potential 
is now monotonic and unbounded from below within the interval £1 < £ < £2- 
If £2 < ^s, where Ig is the proper radius at which V{R{i.s)) = the geometry will be 
regular both inside and out. Let Rg — R{^s)- Rg is given by 

and 

£s = ei+ , ^ . (5.2.2) 

Jr, VI - V{R) 

where V{R) is given by (5.3), (5.4) and Ri is given by Eq.(5.5). The ADM mass, nioo = 
m{l2)- Therefore 

2m^^2Am 87rp2^ (5.2.3) 

-K2 -K2 O 

When I = is, R' = —1 and the QLM vanishes, m{is) = 0. If £2 = is the geometry is flat 
outside. However, before it can become asymptotically flat, the geometry must close in a 
bag of gold. The empty top of the bag consists of a flat cap (in which R' = —1). R{£) = 
at 

is = is + Rs. (5.2.4) 

This is the same (non-generic) kind of singularity we encountered in the uniform density 
model. 

Recall that a regular closed bag requires the integrability condition m{is) = to be 
satisfled. A necessary condition is that p assume a maximum (or a minimum) somewhere. 
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The simplest realization of such a configuration involves three constant density slabs with 
Pi < P2 > Ps (or pi > P2 < Ps)- To satisfy the integrability condition we will be required 
to tune the relative values of these densities and their supports appropriately. The model 
we are considering with pi < p2 with a vacuum exterior represents a special case. The 
value £s given by Eq. (5.2.4) is the unique critical value of £o supporting an everywhere 
regular closed cosmology for the given parameter values. 

If £2 > £s, -R' < — 1 at £ = £2 and m(£) will be negative for all subsequent i until a 
singularity is reached at some point is- If 5 in addition, £2 > is where 

^^ dR 



is = is + 

is finite, the geometry closes with a metric singularity R' ^ —00 before we ever reach £2- 
In the neighborhood of R = 0, R satisfies Eq.(2.8) with singular solution (2.9). The local 
value of the QLM at the singularity is given by m{£s) = Am. Note that the quasilocal 
mass at the singularity is finite and determined completely by the remote behavior of the 
sources. If £2 < ^s-, the interior is fitted with a negative m exterior vacuum cap with the 
same singular behavior. 
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6. APPROXIMATING THE BOUNDARIES IN C 

6.1 Inequalities of Sufficiency and Necessity 

It is clear that a generic energy profile can, in principle, be approximated arbitrarily 
closely by a a piece-wise constant density model along the lines of our examination of the 
two simplest such models in Sects. 4 and 5. Even still, however, it is not obvious how 
to exploit such models to characterize the boundaries partitioning C which were discussed 
in Sect. 3. The approach we will take is to establish necessary and sufficient conditions 
determining the circumstances under which the geometry can possess trapped surfaces 
or singularities [8, 9]. If we could solve the problem exactly, as we did on the piece- 
wise constant density subset, presumably the necessity and sufficiency conditions could be 
formulated as if and only if statements with an appropriate choice of variables.* 

Let us examine in very general terms how closely these necessary and sufficient condi- 
tions identify the (singular / non-singular say) partition S of C. The necessary condition 
will provide us with a surface iSnec bounding some region in the configuration space Cnec in- 
terior to which we are assured that the geometry is non-singular. The sufficiency condition 
will provide us with some other surface iSsuflf bounding a region in the configuration space 
Csuff exterior to which the geometry is always singular. Cnec will be some proper subset of 
Csuff . The set difference Csufi — Cnec will not be empty, representing a grey area containing 
both singular and non-singular geometries. As the constant density model illustrates, this 
discrepancy might be a consequence not only of our inability to solve the problem exactly 
but also of our inability to identify an optimal set of variables. 

Ideally, we would like both the necessary and the sufficient conditions to be sharp, 
namely that a configuration exists saturating the condition, though we do not necessarily 
expect the two to be sharp simultaneously. 

6.2 Sufficiency and Necessity in terms of M and (.q 

Let us rewrite Eq.(2.1) in the form 



* Of course, one could argue that whenever the problem is exactly solvable, in principle 
at least, one could simply tick off geometries with i?' = or i? = at some finite value 
of l. We do not consider this as a solution to the problem because it is telling us nothing 
about the physics conspiring to produce these coincidences. 
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47rpi?2 + (^^/y = ^ (^1 + (i?')^) , (6.2.1) 

and integrate out from £ = up to the surface value Iq: 

M{io)+RR' =r, (6.2.2) 

where we define 



1 



eo 



r = - di{i + {R'y). (6.2.3) 

2 Jo 
We now have two different representations for the first integral of the constraints. Earlier we 
exploited the expression Eq.(2.5)/(2.6) which involves the definition of the QLM. Eq.(6.2.2) 
involves M{£o) and will provide the basis for our derivations of all inequalities involving 

M. 

6.2a Trapped Surfaces: Sufficiency 

In our examination of the constant density star, we found that the sufficient condition 
for the presence of a of trapped surface is as good as one could hope for — the inequality 
(4.11) is sharp. This success is repeated for generic density profiles. 

We begin by eliminating the term involving R'{io) from Eq. (6.2.2) in favor of the 
expansion of the optical scalar G_ = G = Q+ (see appendix of I): = 2R'/R. If the 
surface is not trapped, Q{£o) > and thus M{io) < T. We now require an upper bound 
on r. When p > and the interior is regular, then (-R')^ < 1 so that 

T<io. (6.2.4) 

Thus 

M(4)<4. (6.2.5) 

In other words, if M{£o) > £o the surface must be trapped or the interior must be singular 
in which case it also possesses a trapped surface.* BMO'M have demonstrated that this 



* Note that in the derivation, we did not have to assume that the interior did not contain 
a trapped surface. 
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inequality is sharp. They do this by constructing a model in which M(£o) = ^o(l ~ e) but 
the surface is not trapped. 

It is possible to tighten the sufficiency condition whenever p' < 0. We exploit Eq.(2.5) 
to express F as a functional of m(£): 








\{i~{Rr) 


mm 




R 





(6.2.6) 



We now exploit Eq.(2.10) which holds whenever p' < to obtain 



r<4-|f.«p«^^.„-MM. ,6,2.7) 

which is a tighter bound on T. Thus, in place of Eq. (6.2.5) we have 

^M(4)<4. (6.2.50 

Thus, if M(£o) > 1^0 the surface of the monotonic configuration must be trapped or it 
contains a trapped surface. 

A nice feature of Eq. (6.2.5') is that not only is it sharp, but we also know that 
it is saturated by the constant density star. In addition we note that it is easier to 
form a trapped surface when p is monotonically decreasing. This is to be contrasted 
with the impossibility of finding a singularity when p is monotonically decreasing. If we 
scale the value of p or its support, a critical point will always be reached at which a 
trapped surface forms. However, there is no corresponding point marking the formation 
of a singular geometry. This exposes the inadequacy of the intuitively reasonable notion 
that by increasing p on a constant support we first produce a trapped surface and then a 
singularity. 

We note that the sufficiency condition can be cast in a weaker form in terms of the 
minimum energy density pMin using pMinVb < M(£o) where Vb is the volume of the star. 
In the monotonic star, this will be slightly more useful on account of the fact that pMin 
will be the surface energy density. More useful inequalities involving pMin can be derived 
by exploiting the calculus of variations in the manner of Schoen and Yau [14] . 
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6.2b Trapped Surfaces: Necessity 

In the constant density star, we found tiiat M was a poor variable for formulating 
necessary conditions for the presence of trapped surfaces and of singularities. 

In general, if the surface is trapped, 6(£o) < 0. Then Eq.(6.2.2) implies that M{Iq) > 
r. We now require a lower bound on F. The weak lower bound on F, F > ^, relies only 
on the fact that (-R')^ — which is independent of the positivity of p. We get M{io) > -^ 
on a trapped surface. Thus if 

M(4) < I (6.2.6) 

the surface is not trapped. It is not sharp and it fails to prevent a trapped surface showing 
up in the interior. Compare Eq. (6.2.6) with the marginally better Eq.(4.9) we obtained 
for the constant density model which, in addition, guaranteed that no trapped surfaces 
showed up in the interior. The reason for the latter stronger statement is the fact that we 
were able to exploit the exact solution to feed global information into the surface condition. 
As an extreme example, consider a star with a high density core of radius £i and a large 
low density mantle. We might then have M{li) > £i but M{£o) < IqJI on the surface. In 
other words, the surface information provides no clue as to the interior physics. For this 
reason the condition p' < does not allow us to 'tighten' the necessary condition the way 
it did the sufficiency condition. Technically, the inequality is operating the wrong way. 
What is possible is to tighten the necessary condition in a configuration with a radially 
increasing energy density p' > to 

M(£o)<34/5. (6.2.6') 

This is the opposite extreme to the central high density core. It is harder to produce 
a trapped surface. However, even on this subset, we still cannot provide a strong form 
of necessary condition prohibiting trapped surfaces showing up in the interior. We note 
that the constant density star is also a special case of Eq. (6.2.6'). Eq. (6.2.6') is not sharp 
because Eq.(4.9) is not. We will return to the issue of formulating necessary conditions 
below using pMax instead of M. 

We will find a necessary condition analogous to the constant density result holds when 
we use pMax instead of M. The reason is that pMax taps into global information without 
the distortions introduced by folding K into the definition of M. 
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6.2c Singular Geometries: Inequality of Sufficiency 

Because M does not discriminate between the regular and the weakly singular ge- 
ometries exhibited by the constant density model, we found that instead of a sufficiency 
condition for the presence of singularities, the inequality (4.10) represented a universal 
bound. 

In general M does discriminate between geometries which are strongly singular and 
those which are not. Therefore it can provide a non vacuous sufficiency condition for the 
presence of singularities. This bound was obtained by BMOM [8,9]. On the monotonically 
decreasing density subset, it again provides a universal bound though not a sharp one. 

To reproduce the sufficiency condition, let us suppose that the geometry is not strongly 
singular. Then R'^ < 1, so that T < Iq and Eq. (6.2.2) implies 

M(io)+RR' <io. (6.2.7) 

In addition, R{i) < i everywhere and R' > —1. The surface term is therefore bounded 
from below by —Iq and 

M{io)<2io. (6.2.8) 

If, in particular p' < 0, then this condition represents a universal bound on M in a MSC. 
It is not possible in such a model to raise the value of M indefinitely by increasing p while 
maintaining Iq constant. What happens when we do this is that pRy[.^y^ saturates within 
the star and as a consequence the increase in M{Iq) is at most linear in Iq. We will see 
how this occurs explicitly below in the context of a two density model. In a large outer 
low density region R does not grow with Iq but oscillates between fixed mimimum and 
maximum values. This is folded into the definition of M. 

We can tighten the inequality whenever p' < just as we did the trapped surface suf- 
ficiency condition. Let us do it in a slightly different way which exploits the Schwarzschild 
inequality explicitly. The Schwarzschild radius, -R(£), provides a bound on the quasilo- 
cal mass 2m{£) < R{i) whenever the constraints are satisfied (whether the geometry is 
singular or not). This is a straightforward consequence of Eq.(2.5). 

m-- = -{R) -. 
The right hand side is bounded above by zero {R''^ > 0). Thus 
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m--<0. (6.2.9) 

Equality obtains on an apparent horizon R' = and when R = where 171 = 0. 

The condition (2.8) and the Schwarzschild inequahty (6.2.9) together then imply that 
pR^ < 3/(87r) inside matter. Thus 

M{io) =47r f " diR^p< ^£o- 
Jo ^ 

It is curious that, unlike the corresponding trapped surface condition, this inequality is 
not sharp. 

6. 2d Singular Geometries: Necessity 

Suppose that the geometry possesses a singular surface at Iq. Then R{£q) = 0. What's 
more the product RR' appearing in Eq. (6.2.2) also vanishes at a singularity even though 
R' itself might diverge. Thus 

M{io) = r > y (6.2.9) 

which is the same as the necessary condition for the formation of trapped surfaces. This 
condition is weak — it does not even demand the positivity of p. Neither does it say 
anything about the interior. 

The necessary conditions we have obtained are far weaker statements than we would 
like, failing to eliminate the possibility of trapped surfaces or singularities lurking in the 
interior. What is clear by now is that M{io) and io are not the ideal variables for this 
purpose. 

We turn now to the task of formulating alternative necessary conditions. 

6.3 Necessity in terms of pMax and io 

In the context of the constant density star we discovered that the inequalities which de- 
scribe the trapped/untrapped, singular/non-singular partitions of geometries is expressed 
with greater precision in terms of p and £, the more so for the inequalities of necessity. 
One reason why this is so is that the single constant value encodes global information. 
However, another reason is that the inequalities do not involve R. What therefore are the 
appropriate generalizations of the constant p which are independent of R7 
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The generalization which permits us to formulate necessary conditions is the maximum 
value of p on the support of matter, pMax- The dimensionless combination of pMax and Iq 
is PMax^o- Such conditions can be cast in the form: trapped/singular on surface or in the 
interior implies pMax^o — some constant. What's more the constants will be different in 
the two cases just as they were for a constant density star. Not only are these necessary 
conditions stronger than those obtained earlier in terms of M for characterizing the surface, 
they also detect trapped surfaces/singularities in the interior if they are present. 

6.3a Singular Geometries 

If the interior geometry is singular then R{ii) = at some io > £i > 0. We integrate 
Eq.(6.2.1) from £ = up to £i. The integral over the divergence results in a surface term 
which vanishes due to the boundary condition on R ai £i. We obtain 



/ M{1 + R''^) =8n dipR^ . (6.3.1) 

Jo Jo 



It is crucial for the successful implementation of the integration that the product RR' 
vanishes at a singularity even though R' itself might diverge. We now apply the Holder 
inequality with the L°° norm on p to obtain the bound 

£i I'll 

dipR^ < pMax / diR^ (6.3.2) 

'0 Jo 

on the left-hand side. The bound is clearly exact if p is constant but is a weak estimate if 

p varies a lot. Hence we get 

I dl{l + R'^) < SyrpMax / diR^ . (6.3.3) 

io io 

The nice thing about this expression is that we can exploit a Sobolev inequality to get a 

second bound on the quantities occuring in it. For functions on the finite interval (0,£i) 

which vanish at the endpoints, there is a positive constant 5*0 such that 

So / diR"^ < / diR'^ . (6.3.4) 

Jo Jo 

The inequality is saturated by the trigonometric function R{£) = sm{n£/£i) which also 

determines the optimal value of 5*0: 5*0 = t^^/£i- Sq is just the ground state energy of 

a quantum mechanical particle in a box of width, ii. In a constant density star R{£) is 
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exactly of this form so that the second estimate involved here is also exact. Substituting 
Eq.(6.3.4) in (6.3.3) we get 

rh ^2 nil 

/ M = h< (SyrpMax - yr) / ^^^' ■ (6-3.5) 

While we do not have a lower bound for R' , we do have an upper bound which is that 
R' <1. This then gives us that R < i and that 

diR'^ < -^ . (6.3.6) 

3 

When this is substituted into Eq. (6.3.5) we get 

PM.Ji > ^ + |- • (6.3.7) 

8 on 

Hence we conclude that if pMax^i < tt/S + 3/(87r) we will not reach a singularity within a 

proper distance ii from the origin. This estimate is only one of a family of such estimates 

that can be derived. Let us begin by multiplying Eq. (6.2.1) across by i?" where a is some, 

as yet undetermined, constant. We can write (6.2.1) as 



R" + (2a + l)i?"(i?')^ - 2{R''+^Ry = 87rpi?"+2 . (6.3.8) 

We integrate this out to a singularity at £ at £i. The divergence gives us a surface term 
which will vanish so long as a > — ^. Thus we get 



Jo 



"h , , r^i 



f ' diR- + i^^^±^ / ' dlilR'^yf = 87r / ' rf£pi?"+2 . (6.3.9) 

(a + 2)^ Jo Jo 

We now repeat the argument as above, replace p with pMax and use the Sobolev inequality 
(with the same constant) to get 






^iJo 



■^^l'^2 ^ TT , ,. p„+2 



di{[R—yy>^l diR^^^ (6.3.10) 



When these are substituted into Eq. (6.3.9) we get 



■i 
dlR" < 





87rpMa 



4(2a+l)7r 



(a + 2)2£2 
The last estimate needed is a bound on the ratio 
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R"^^ . (6.3.11) 



rf£i?"+7 / MR". (6.3.12) 

'Jo 

The very crude bound, l\, follows immediately from the fact that R' < 1 which gives 

us that R < £i. Because R appears in the denominator, we cannot naively exploit the 

bound R < £ as we could in Eq. (6.3.6). One can, however, show that the linear function 

R = i maximizes the ratio (6.3.12) for all functions which satisfy the boundary condition 

-R(O) = 0, and the two constraints R' < 1 and -R > (see appendix II). Thus 

^ < 77^^? , 6.3.13 

J^'diR- 3 + ai' 

holds. We thus obtain 

,2 ^ (l + 2a)7r 3 + a 1 /^ o . .x 

The dominant term on the RHS is the first. The best value for it occurs when a — 1. We 
get 

PMaJi >l + ^- (6.3.15) 

6 47r 

The RHS of Eq. (6.3.15) is slightly larger than that of Eq.(6.3.7). Note that a = 1 cor- 
responds to an integrand RR'^ which tends exactly to a constant at a singularity. More 
accurately, the RHS of (6.3.14) is maximizes when a ~ 0.8. This is less than 1% better 
than Eq.(6.3.15). None of these numbers are particularly impressive, they are, at best, 
only half the value we get in the constant density star. 

Despite our effort, it is clear that the estimate Eq.(6.3.13) is the weakest link in the 
argument. 

6.3b Apparent horizons 

Suppose now that the geometry possesses an apparent horizon, R'{£i) = at some 
^0 > ^1 > 0. We again integrate Eq.(5.2) from £ = up to £i: As before, the boundary 
condition kills the surface term and we are left with an identical equation: 

de{l + R'^) = Stt / dlpR^ . (6.3.16) 

Jo 
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The first approximation proceeds identically 

de{l + R'^) < SyrpMax / MR^ . (6.3.17) 

io 

For the second approximation, we again exploit a Sobolev inequality involving an appro- 
priate function space on the right hand side. If R{t) is a function on the finite interval 
(0, £i) such that -R(O) = and R'{£i) = 0, then there is some constant Si such that 

rii rii 

Si / diR^di< / diR'^. (6.3.18) 

Jo Jo 

2 

The inequality is again saturated by a sine function, R{i) = sin(7r£/2£i) giving Si = j^. 
In the case of a constant density star, both approximations are exact. So 

/ di = ii< (STrpMax - j^) / diR^ . (6.3.19) 

Jo ^^i Jo 

We again use Eq. (6.3.6), and we finally get 

PMax^? > ^ + ^ • (6-3-20) 

This is a much better estimate than the corresponding bound (6.3.7) or its improvements 
for the singularity, largely due to the fact that (6.3.6) is much sharper in this case. In fact 
it is about 3/4 of the number that we get in the constant density star. 

An expression equivalent to Eq.(6.3.14) can also be written down. We have that 

, .2 ^ (1 + 2a) TT 1 3 + a 

pMax«i > — — — ^- + — -^. 6.3.21 

[2 + ay 8 87r 1 + a 
The maximum of this expression occurs at about a = 0.25 and this improves on the bound 
Eq. (6.3.20) again by about 1%. 

If p is bounded from below within the star, its minimum value, pMin can be exploited 
to provide alternative sufficiency conditions. However, if p vanishes anywhere this fails. 
In addition, because the inequalities of sufficiency are derived by asssuming that the star 
does not possess a trapped surface/singularity we cannot exploit Sobolev inequalities in 
the way did for the necessary conditions. 
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7. A LOWER BOUND ON THE BINDING ENERGY OF A SPHERICALLY 
SYMMETRIC DISTRIBUTION OF MATTER 

7.1 Momentarily Static Spherically Symmetric Shell 

An important special case of the two-density model is the shell with* 

p{(.)=a5{(.-e.o). (7.1) 

This model will be useful in two regards. Firstly, it represents the configuration with least 
binding energy. Secondly, when we generalize our analysis to configurations with non- 
vanishing current, it provides a useful exactly solvable model in which we can examine the 
behavior of the geometry near singular points [3] . 

Inside the shell, space is flat so that R = L As a consequence the material energy M 
is equal to its Newtonian value 

M = 4nail . (7.2) 

By integrating Eq.(6.2) across £ = £o we flnd that R' suffers a discontinuity at the shell 
given by 

AR' = -Anaio . (7.3) 

Outside i?'(£o+) = 1 + AR' . If a is positive i?'(£o+) will be bounded above by one. As we 
have seen M does not see the discontinuity in R' . The ADM mass nioo, however, does. 



nir 



m{eo+) = j{l-R'\£o+)). (7.4) 



A horizon must form at some point outside the shell if 4na£o > 1. 
If 4Tva£o < 2, R' > —1, nioo is positive and the geometry is regular. 

When 47rcr£o = 2, R' = —1, rriinfty = and the geometry is weakly singular. At this 
value, i?'(£o+) = —1, the ADM mass m vanishes and the geometry outside is flat and 
closed as well. The bag of gold consistes of two flat caps sewn together along the shell. 



* In the notation of Sect. 5, the shell corresponds to pi = 0, and the limiting value 

o- = P2{h -£i)- 
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If Aivalo > 2, R'{£qj^) < —1, rrioo is negativeand the geometry suffers from a strong 
singularity in the same way as it does in the two-density model discussed in Sect. 5. 

A useful expression for rrioo is obtained by writing AR' in terms of M and Iq. We use 
(7.2) to express (7.3) AR' = —M/£q so that Eq.(7.4) can be rewritten in the form 

m = M---. (7.5) 

Zto 

The binding energy M — m = M^/(2£o) is the Newtonian value. The distribution of matter 
with the least binding energy in Newtonian gravity is a shell. 

7.2 Lower Bound on the Binding Energy 

It was conjectured by Arnowitt, Deser and Misner but only proven recently by BMOM 
that, in general, [15, 10] 

M-m>—. (7.6) 

This provides a lower bound on the binding energy of a spherically symmetric distribution 
of matter of fixed M and Iq. The inequality is sharp. The distribution of matter which 
saturates this inequality is the shell (compare Eq.(3.17)). 

Note that if I is decreased while M is kept fixed, the binding energy increases — the 
more compact the material system the larger its binding energy. BMOM's proof relied on 
the use of a conformally fiat coordinate system. We will reproduce the proof using the 
proper radius directly. The trick, as demonstrated by BMOM is to mimic the calculation 
for the shell. We rewrite Eq. (6.2.2): 

We now substitute for i?' in the expression (2.5) for m: 

m = |(l-i?'^)=-^ + |(l-F2) + F)M, (7.8) 

where we introduce the dimensionless quantity F := -^. We now express the difference as 
the sum of the lower bound we wish to establish and a remainder, Q 

M-m= — + Q, (7.9) 
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where 

A/f ^ / 1 1 \ T? 

Q = — (^ - ^) + (1 - i^)M + |(F2 - 1) . (7.10) 

We will show that Q is always positive. The key to doing this is to note that if the geometry 
is non-singular, then 

1<F<|. (7.11) 

The upper bound is simply Eq.(6.2.4).* To obtain the appropriate lower bound, let us first 
examine the difference between the proper radius and the circumferential radius 



i-R= M{l-R'). 
Jo 

We can place a lower bound on this difference using the inequality 1 + R' < 2: 

1 r^' 

e-R>- / di{i-{R'f). 

2 Jo 
This can be inverted to provide the lower bound 

- / di{R'f>R--. 

From this we can deduce that F > 1. We therefore have established both an upper bound 
and a lower bound on F. 

We can now demonstrate that Q is always positive. We note that both the first and 
third terms are positive. The discriminant of Q is given by {F — 1)IF + 1 — 2i/Rj and is 
negative. The quadratic therefore possesses no real root. Thus Q is positive everywhere. 

Note that Eq.(7.6) implies m > iS M < 2£. This reproduces the sufficiency condi- 
tion, Eq.(6.2.8). 



* We also have already got a weak lower bound on F which does not rely on the positivity 
of p, i/{2R) < F. 
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8. CONCLUSIONS 

In this paper, we have attempted to identify generic features of asymptotically flat 
spherically symmetric solutions to the Hamiltonian constraint when the spatial geometry 
is momentarily static. 

Our focus has been on the characterization of the two strong field features of such 
geometries, apparent horizons and singularities. The simple exactly solvable models con- 
sisting of piecewise constant energy density profiles provide a useful guideline for the choice 
of appropriate variables to characterize these features. Sufficiency conditions describing 
their existence can be cast as inequalities between M and Iq. Matching inequalities of ne- 
cessity are better cast in terms of pMax and £o- We found that the latter inequalities could 
be improved using simple functional inequalities. Because these inequalities are not sharp, 
it is likely that they can be improved with a more judicious exploitation of these inequal- 
ities. This work might provide clues towards the identification of appropriate variables 
with respect to which to formulate the necessary part of the hoop conjecture [11]. 

The techniques introduced in this paper prepare the ground for the examination of 
the constraints when Kab 7^ 0. In paper lib, we will examine the analogous problem in 
this more general case. 

At the end of paper I, we introduced the optical scalar plane as a representation 
of the phase space of the theory. In this representation non-singular solutions of the 
constraints appear as bounded closed trajectories each containing the point P = (2, 2). The 
moment of time symmetry solutions we have considered in detail in this paper correspond 
to degenerate trajectories that run along the diagonal, a;_|_ = uj- starting out at the point 
P. If the solution possesses a trapped surface it will cross the origin, (0, 0), corresponding 
to R' = at least once. If it is non-singular, it will make an even number of crossings 
before returning to the point P. In particular, we saw that when the energy density profile 
is monotonically decreasing, a configuration can cross the origin an arbitrarily large even 
number of times. If the solution is singular it will cross the origin an odd number of 
times before ultimately proceeding towards the point, Q = (—2,-2). At this point, the 
QLM vanishes and R' = —1. Only in situations of zero measure occuring when the energy 
density is fine-tuned in such a way that the QLM tends monotonically to zero will the 
trajectory terminate at this point. In general, once Q is breached the trajectory cannot 
re-enter the bounded interval along the diagonal between Q and P. What must occur is 
that it continues monotonically along the diagonal towards unboundedly large negative 
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values of R' . 

When Kab 7^ 0, there are many more possibihties [3]. Solutions may possess a future or 
a past apparent horizon. These horizons will no longer generally coincide with the extremal 
surfaces of the three-geometry. Neither does the existence of the former necessarily imply 
the existence of the latter. 

Whereas in the MSCs, the only approach to singularity is through the point Q, there 
is now a wide range of possibilities. In general, a non-vanishing extrinsic curvature can 
lead to more severe singularities than those encountered at a moment of time symmetry. In 
addition, the constraints no longer imply that the scalar curvature is finite. The converse 
of the positive QLM theorem is no longer true. Indeed, m can be positive everywhere and 
yet the geometry be singular, m will, however, always be finite if the sources are finite. 

What is remarkable is that it is possible to generalize the necessary and sufficient 
conditions examined in this paper with the appropriate generalizations of M and pMax 
and an appropriate gauge. Not surprisingly, the gauges that do work are procisely the a - 
parametrized linear extrinsic curvature gauges introduced in paper I, with a in the range 
0.5 < a < 00 [1]. 

The bags of gold which occur behind singularities are physically disconnected from 
the exterior and, as such, of little more than curiousity value in the present context, at 
least, at the level of the classical theory. In cosmology, however, the bag of gold can be 
interpreted as a closed universe. 

If £ = is the north pole, £ = £5 is just the south pole of the closed universe. This 
closed universe will generally not be regular. In fact, the argument presented in Sect. 2 
suggests that regular closed universes constitute a subset of zero measure in the set of all 
closed universes. This universe will be singularity free if and only if the QLM vanishing at 
one pole also vanishes at the other pole. We can think of the value M{ls) assumed by M 
at the south pole as a measure of the binding energy of this regular closed universe. 

Using Eq.(2.6'), the integrability condition m{is) = on a regular bag of gold implies 

diR^p' = 0. (8.1) 



A regular spherically symmetric closed universe is impossible if p is strictly monotonic. 

p must possess a maximum or a minimum away from the poles. We will examine the 

configuration space of regular spherically symmetric cosmologies in a subsequent paper 

[16]. 
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APPENDIX I: THE RICCI TENSOR FROM THE SCALAR CURVATURE 

It is always possible to construct the Ricci tensor once we know the scalar curvature. 
We will do this for a general spherically symmetric geometry by appealing to a mini- 
superspace Hilbert variational principle: Let 



d^x^TZ . 



(a.l) 



Then 
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ah 



(a.2) 



^/g Sgah ' ' 2 

To implement the variational argument we need to restore the metric coefficient C^ (defined 
in I) so that the spatial coordinate system is only defined up to spherical symmetry. This 
is because we require as many independent metric components as there are independent 
components of TZah which is two in this case. Up to a divergence which is irrelevant for 
our purposes (here ' refers to differentiation with respect to r), 



Thus 



We now note that 
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Having taken all variations, we can safely set £ = 1. We use Eqs.(a.2), (a.4) and (a.5) to 
express the non-vanishing components of the Ricci tensor in the form 

n-^\n + ^,{R--i)- n'' = ^n + ^. (a.e) 

We now eliminate the second derivative in the second term in favor of TZ and lower deriva- 
tives of R. Eq.(2.12a & b) for the scalars TZc and TZr are then given by the appropriate 
projections. 

In any regular solution of the constraints TZ^ > 7^/4, and TZc < TZ/2 everywhere. 

We note that there is no non-trivial regular everywhere solution of the Hamiltonian 
constraint in which either TZfi or TZc vanishes. In the case of TZ^ = this is obvious 
on account of the above bound. In the latter case, consistency with the constraint (at a 
moment of time symmetry) then implies that 

1 r^ 

pR^ = - pR^R'di , 



-R JO 
so that p ~ R~^ which is singular. 

APPENDIX II: PROOF OF THE BOUND (6.3.13) 

In this appendix, we prove the validity of the bound, Eq.(6.3.13) in the text for all 
a > — 1. For simplicity, we will first provide a proof for a = 1. The generalization will 
then be clear. We define 

F{RJi)= R^M Rdi. (6.1) 

We will prove that 

F{R,h)<^l, (6.2) 

is true, for functions that satisfy 

R{0) = 0, R'{i)<l R>0. (6.3) 

Let us calculate the first variation of F 
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6F^ {3R^ - F)5Rdi / { RMf. (6.4) 

Jo ' Jo 

We claim that the maximum is achieved hy R = L 

In this variational problem, we can only consider variations which satisfy Eq.(b.3), 
in other words we start with a function Rq which satisfies Eq.(b.3) and look for a family 
of functions i?t, such that Rt=o = Ro and that for some range of t, ^ At, R{t) also 
satisfies Eq.(b.3). The 5R in Eq.(b.4) must be the first derivative of such a sequence Rt 
with respect to t, evaluated at t = 0. 

In particular, the allowed variations around R = i must satisfy 

SR' < 0, SR{0) = , (6.5) 

which forces 6R to be negative and monotonically decreasing on the interval (0, L). Now 
let us consider the integral in Eq.(b.4), with R = i and F = ii /2 if we replace 5R with a 
negative constant —C. It is clearly negative, equalling —C£i /2. Notice that 3R^ — F is 
negative in the range (0, ii/V^i) and positive on the rest. This means that, since 5R must 
be negative, the integrand in Eq.(b.4) is positive on the interval (0, £i/\/6) and negative 
elsewhere. 

Since SR is monotonically decreasing if we replace 5R in the integral by its value at 
i = i\l \f^ we make the integral more positive. However, it is still negative. Therefore the 
actual integral is negative and all the allowed variations reduce F . Hence i? = £ is a local 
maximum. 

This is not enough, however. We need to show that no other maximum occurs. To do 
this we note that F < -Rmaxi where -Rmax is the maximum achieved by R on the interval 
(0,£i). This can be shown by extracting -R^^x from the R^ integral. We immediately get 
that {3R^ — F) must be positive near the maximum of R. If we further have that R' < 1 
at -Rmax we can find a positive variation dR, localized near the maximum of R, which 
satisfies the R' < 1 condition. Such a variation will increase F. Thus, in particular, R 
cannot possess an interior maximum. The maximum must be assumed at i = £i and at 
this maximum, R' = 1. 

Thus the only case we have to consider is the case where R achieves its maximum at 
ii and R' = 1 there. We need only to show that if -Rmax < ^i, there exists a variation 
increasing F. 
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Suppose that -Rmax < ^1- We cannot have i?' = 1 on the whole interval, for this would 
imply R = L Therefore we have a point £2 where R{£2)' < 1- An allowed variation is then 
SR = on the interval (0,^2) and SR = 1 on (£2,^1)- For £2 close to £1, the integral in 
Eq.(b.4) can be approximated 

[3{i -ii+ i?„ax)' - F]di =i?Lx - (^2 -il+ i^max)' " F{ii - I2) 

= (4 - ^2) [(4 - ^2? - 3i?„ax(^l - ^2) + 3i?Lx - F] ■ 

(6.6) 
We know that we can choose £2 such that ^1—^2 < -Rmax because of the positivity of 
R. But this, in conjunction with F < -R^^x' i^ sufficient to show that the expression in 
Eq.(b.6) is positive. This completes the proof. 

The generization of this proof is straightforward. We can prove that 

^„,«,,,.i^<(i±i!)V. (.7) 

/o ' R'^M 3 + a 

for any a > — 1. The technique is as before, first of all we prove that R = I is a, local 
maximum, then we prove that if -Rmax < (-1 there exists an allowed variation which increases 
F. 

We have, 

fn V[2 + a]R^+'' - aFR-^+'')6Rde 
SFa = ^^^-^ . . 6.8 

If a < we immediately see that increasing R increases F. Therefore R = i is the 
maximum. We can now restrict our attention to the case where a > 0. 

If we replace dR in Eq.(b.8) with a negative constant —C we can do the integration 
and show that it is clearly negative, equalling —2C£i /(3 + a). It is clear that the allowed 
6R^s must be negative and monotonically decreasing. The integrand in Eq.(b.8) is positive 
in the range (0,^2) and negative in (£2,^1), where £2 is given by 

2_^aa + aL 2 . 2 
^2-(2 + a)(3 + a)^^ < ^^ " 

Therefore we can replace 5R by its value at £2 and make the integral more positive. 

However, since the integral with a negative constant is negative, the correct integral is 

even more negative. Therefore R = £ is a local maximum. 
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The rest of the proof is just as before. We know that F < -R^^x ^o therefore the 
quantity {[2 + a]R^~^"' — aF R~'^~^°') is positive near the maximum of i?. Therefore a positive 
variation, locahzed near the maximum, increases F. The only situation we need to deal 
with is where the maximum of R occurs at i = £i and simultaneously R' = 1 there. Let us 
find the point £ = h where R'{£2) < 1 and R = £ — £i + -Rmax in the interval (£2, £i)- Now 
an allowed variation is one which is zero on the interval (0, £2) and a positive constant on 
the interval (£2,^1)- For such a variation we can do the integration in Eq.(b.8) to give 



6Fa{ / R^'Mf =R^+,l - {£2 + i?n.ax - 4)'+" " Fi^^ax + ^(^2 + i^max " £lY ,, ^, 

Jo (^-9) 

=^max[^max - F] - {£2 + i?niax - ^l)"[(^2 + ^max - £l) - F] . 

The first term in on the second line Eq.(b.9) is positive and also the first term is greater 
than the second term because -Rmax > £2 + -Rmax — £i- Therefore the variation of F is 
positive. 
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Figure Captions 

Fig. (4.1) M{io)/iQ vs. io Whereas M is monotonic in io by construction, tlie ratio M/ig 
is not — even before tlie geometry becomes singular it begins to oscillate. The ratio rises 
monotonically for small values of io, regardless of the formation of an apparent horizon 
until a maximum value is reached when the surface becomes a decreasing point of inflection 
of R{i) ■ Contrary to what we might have expected, the maximum does not obtain with the 
appearance of a singularity. As the star is made even larger, the ratio begins to decrease 
until it re-assumes the value 3/4 as the surface becomes singular. Subsequently, the ratio 
oscillates with decreasing amplitude with one complete oscillation between each occur ance 
of a singularity. As £o ^ oo, the ratio tends to the constant value 3/4 independent of p. 
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